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Abstract. 

We consider the optimal covering of fractal sets in a two-dimensional space using 
ellipses which become increasingly anisotropic as their size is reduced. If the semi- 
minor axis is e and the semi-major axis is 5, we set S = e a , where < a < 1 is an 
exponent characterising the anisotropy of the covers. For point set fractals, in most 
cases we find that the number of points N which can be covered by an ellipse centred 
on any given point has expectation value (M) ~ e" , where f3 is a generalised dimension. 
We investigate the function /3(a) numerically for various sets, showing that it may be 
^ different for sets which have the same fractal dimension. 

00 
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1. Introduction 

O 

Figures CD and [2] illustrate two different fractal point sets in two dimensions (their origins 
will be described shortly). The fractal dimensions (as defined in [IJ[2]) of these sets are 
very similar (the correlation dimension are ~ 1.76 and D2 ~ 1.71 respectively), 
but inspection of these figures suggests that the fine-scale structure of these sets is 
quite different. In this paper we consider whether it is possible to make a distinction 
between the local structures of fractals by using covering sets which become increasingly 
anisotropic as the covers are made smaller. We show that this can distinguish different 
sets which have the same fractal dimension. Moreover, we shall argue that the use of 
anisotropic covers can be important in analysing the scattering of radiation from fractal 
distributions of matter. 

The sets illustrated in figures [1] and [2] are point set fractals which arise as models 
for the distribution of particles resulting from two distinct physical processes. Figure [1] 
illustrates a model of particles in a turbulent flow, in circumstances where inertial effects 
are large enough to ensure that the particles are not simply advected. It is known that 
such systems may exhibit clustering [3] and that the attractor is a fractal measure 
[U |5]. An enlargement of a subset of figure [1] shows that the particle distribution is 
locally highly anisotropic, with the particles tending to be clustered close to lines in the 
plane, even though the statistics of the distribution are globally rotationally invariant. 
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Figure 1. Distribution of particles with significant inertia moving in a two-dimensional 
area- preserving velocity field (the equations of motion are specified in section |3J) . The 
particles tend to cluster and for the parameters used in this simulation (see section \S§ 
the correlation dimension is D2 ~ 1.76 



This is consistent with the structure of strange attractors in low- dimensional dissipative 
dynamical systems, where the attractor has a local structure which is the Cartesian 
product of a line and a one-dimensional Cantor set [6]. The second example, figure El is 
a diffusion-limited aggregation (DLA) cluster in two dimensions |7] . The local structure 
of this set cannot be approximated as a Cartesian product. Also, if a subset of a DLA 
cluster such as figure El is examined at a higher magnification, it does not appear to be 
more markedly anisotropic in structure. Thus it appears as if some fractal sets, such as 
that illustrated in figure (H may exhibit a pronounced anisotropy under magnification, 
whereas this effect is much less pronounced (or possibly absent) in others, such as figure 

El 

Throughout this paper we confine our attention to point sets in the plane, but 
generalisations to higher dimensions are obvious. We use the following approach to 
characterise the local structure of a point set. Take a given element of the set, and 
consider an ellipse centred on this point, with its semi- minor axis of length e and its 
semi-major axis of length 5 = e a , where < a < 1. We assume 5 <C £ and e <C £, where 
£ is the characteristic lengthscale of the system below which fractal scaling is observable, 
and that e ^> eo, where eo is the lengthscale where the fractal scaling is cut off by the 
finite number of points sampling the fractal measure. We then choose the orientation 
of the cover so that it maximises the number of other points which are contained in 
this ellipse. We denote the number of points under this optimally-oriented cover by M . 
We repeat this for ellipses centred on other randomly selected points in the set, and 
compute the average value (N) of the number of points which can be covered. In most 
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Figure 2. Distribution of particles determined by a diffusion-limited aggregation 
(DLA) process. This distribution has (approximately) the same fractal dimension as 
figure [0 D 2 » 1.71 



of the examples of point-set fractals which we investigated, the mean number of points 
in this ellipse is found to have a power-law dependence: 

(M(e, a)) ~ , 5 = e a (1) 

where the exponent (3 depends upon a. In the case where a = 1, the ellipse is a circle, so 
that this case reduces to a definition of the correlation dimension: D 2 = (3(1) (compare 
with the discussion of the correlation dimension in [8]). In general f3 must decrease 
monotonically as a decreases. 

A related approach to the characterisation of fractal sets was proposed by 
Grassberger and co-workers [H [H [10l [11], who considered covering a set (which is 
embedded in a d- dimensional space) with <i-dimensional ellipsoids, with principal axes 
€i, ordered so that e\ > €2 > ■ ■ ■ > e^. Their work emphasises the case where the local 
structure of the fractal is a Cartesian product of sets, with dimensions Aj, ordered so 
that Ai > A2 > . . . > Ad- It is asserted that the ellipsoids cover most efficiently when 
they align with the principal axes, such that the longest axes align with the directions 
of the highest density sets. According to this hypothesis the number of points covered 
is expected to satisfy 

(AO-eNf-ef 3 -.^. (2) 

Examining the dependence of Af upon the e« would allow the partial dimensions Aj 
to be determined. This approach was mentioned in several papers [8], [9J [101 ED] > with 
the motivation to characterise a fractal set by means of its partial dimensions, Aj, 
satisfying Y^i=\^% = D 2 . These works do not prescribe how (or whether) the ratio 
6j+i/ej approaches zero as e\ — > 0. In our work this limiting behaviour is specified by 
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Figure 3. Scattering of a beam of light with wavevector k, width 5 and wavelength 
e: a cluster of J\f particles covered by an e-6 ellipse scatters light coherently when the 
condition for specular reflection is satisfied. The wavevector k' satisfies the condition 
for specular reflection if the major axis of the ellipse is perpendicular to k — k' . In this 
case the scattered intensity from the particles under the cover is increased by a factor 
of order N. 



the parameter a. Our numerical investigations encompass fractal sets which have a 
Cartesian product structure, and those which do not. 

In the case of point set fractals which represent a physical distribution of particles, 
the existence of a highly anisotropic local structure has important physical implications 
for the scattering of light. It is known that s-wave scattering of light from a fractal 
point set gives rise to an algebraic relation between the scattering wavenumber k and 
the scattered intensity /, such that 

I(k) ~ k~° 2 . (3) 

Heuristic arguments supporting this are developed in [12]; see [13] for a discussion of 
some of the complications which can arise in verifying this suggestion. The distribution 
of this scattering intensity may, however, be highly inhomogeneous. In particular, if 
the particles have a strong tendency to accumulate along lines in two dimensions (like 
the example shown in figure [1]), or on planes in higher dimensions, light may scatter 
specularly from these structures. We argue that this motivates the investigation of the 
anisotropic covering with ellipsoids. Consider weak scattering of light with wavelength 
e which propagates as a beam of width 5. When the path length for light scattered 
from different particles is large compared to e, the scattering of light from M particles 
is incoherent, so that the contribution to the scattered intensity is I ~ N '. If, however, 
an ellipsoid of major axis S and minor axis e can be aligned to cover TV particles, then 
there will exist directions where the path length difference is less than one wavelength, 
so that this set of /V particles scatters light coherently (see figure [3]). In these directions 
where the condition for specular reflection is satisfied by the optimal covering ellipsoid, 
there is a greatly increased intensity / ~ TV 2 . 

We have argued that light scattering from fractal sets may be extremely 
inhomogeneous because of specular scattering from particles which accumulate on 
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Figure 4. In the case where the fractal set is a Cartesian product of two Cantor sets, 
with dimensions D x and D y in the x and y directions, the optimal covering ellipse 
might be expected to have its major axis aligned with the direction corresponding to 
the denser set. In the case D y > D x , the major axis of the ellipse would be expected 
to align with the y-axis (vertical in this figure). Here we illustrate a sample of the 
actual optimal covers, which have a distribution of angles of their principal axes. 

surfaces. This motivated us to study the covering of a fractal by anisotropic covering 
sets, and the results are reported here. We investigate the dependence of the generalised 
dimension /3 upon the anisotropy exponent a. We show that the form of the function 
(3(a) can distinguish between different fractal sets which have the same value of the 
correlation dimension D2 = The description of light scattering proved to be a 

very complex issue, which will be considered elsewhere. 

2. Some elementary estimates 

Before discussing our numerical investigations, we consider some simple arguments 
about the form of the generalised dimension (3(a), defined by equation (JTJ. 

First we address the issue of whether the exponent exists. When a — 1, the 
exponent (3 coincides with the correlation dimension of the set. For other values of a, 
we do not know of any general argument which proves that the dependence of the optimal 
covering M has a power law relation to the size e of the covering elements. For most 
of the point-set fractals which we examined we did find good numerical evidence that 
(M) ~ e 13 for small values of e (extending down towards values of e where the discrete 
sampling of the set becomes a limitation). The exceptions occurred for some cases of 
the motion of inertial particles in a random flow, and for some of the sets considered in 
section 0], including the example illustrated in figure HI these are discussed in sections [3] 
and H] respectively. 

Next we give an upper bound on (3(a). Let us consider the expectation (Af) for 
the case where M is independent of the orientation of the ellipse. The total number of 
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particles in a disc of radius 5 is proportional to 5° 2 . If the set is locally isotropic then 
the coverage is independent of the angle of the major axis of the ellipse, and a fraction 
e/S of these points lie in the ellipse. Recalling that 5 = e a , for this isotropic fractal we 
have (A/") ~ S^'h ~ e 1 - Q+aD2 = e^ a \ so that 



is an upper bound on the exponent f3. 

We were not able to obtain a precise and non-trivial lower bound for /3(a), but 
the following argument suggests a possible form for a lower bound. Consider the case 
where the point fractal samples a Cartesian product of two one- dimensional Cantor sets, 
with dimensions D x and D y . In the following we assume that D y > D x . Because the 
set is 'denser' in the direction of the 7/-axis, we expect that the optimal alignment of 
each ellipse is when its major axis is aligned with the y-axis (the same hypothesis was 
proposed in [EJ |9]). The expected number of particles captured by a covering ellipse is 
then (J\f) ~ S Dy e Dx = e aD y +D * ( so that the dimension of this product set is 



Now consider the smallest possible dimension which could be achieved according to this 
argument, if we allow the dimensions D x , D y of the component sets to vary so that 
D x + D y = D2. The greatest number of particles in the ellipse, and hence the smallest 
dimension, is obtained by setting D y — 1 and D x = D 2 — 1. This gives a putative lower 
bound to the dimension: 



The assumption that the covering ellipses align precisely with the y-axis is not really 
correct, as evidenced by figure HI It is plausible that the probability for an ellipses to be 
significantly mis-aligned decreases as e — > 0, but we were not able to obtain conclusive 
numerical evidence. We note that the argument leading to the proposed lower bound, 
/3_(a), is very similar to that presented in [8j [9] to motivate the concept of partial 
dimensions. 

In addition to the fact that the optimal covering ellipses do not align perfectly 
with the preferred axes, there is a further complication which could affect the argument 
leading to the estimate in equation ([6]). In sets such as that illustrated in figure HJ the 
local structure is a Cartesian product of a line and a one-dimensional Cantor set. The 
line is, however, curved. It is interesting to consider whether this curvature can alter 
the estimate in equation ([6]). In order to consider the effect of this curvature, introduce 
two coordinates: y is a coordinate for the expanding (unstable) manifold centred on the 
reference particle and X is a coordinate for the stable manifold. Because the manifolds 
are curved, the equation of the unstable manifold will be of the form x ~ Cy 2 for small 
y, where C is a constant. Consider a family of ellipses, with fixed S and decreasing a 
which are aligned so as to provide the optimal cover for a cluster centred on a reference 
particle. As a is reduced, the curvature of the unstable manifold can become important, 
because it will take a cluster of particles outside the covering ellipse (see figure EJ). This 
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Figure 5. In many dynamical systems the attractor is locally a Cartesian product of 
a Cantor set and a line. If this line is curved, an ellipse may cease to be an efficient 
cover as a is decreased. The x- and y-coordinates are referred to in section [2l 

happens when x(S) ~ e, that is when Ce 2a ~ e. Accordingly, we might propose that the 
optimal covering strategy of aligning the principal axis of the ellipses with the unstable 
manifold starts to break down when at a critical value of the exponent, a c = |. On 
the basis of this argument we would expect that (3(a) might exceed equation (jSJ) when 
en < ~, in cases where the fractal is generated by a dynamical system for which the lines 
representing the unstable manifold are curved. 

3. Numerical investigations of dynamical fractals 

Figure [1] illustrates the distribution of independently moving inertial particles in a 
random velocity field. The equations of motion for the position of a given particle 
are [3] 

r = v, v = — 7[« — u(r(t), t)] (7) 

where 7 is the rate at which the particles relax towards the fluid velocity, and where 
u(r,t) is a randomly fluctuating velocity field satisfying the incompressibility condition 
V • u = 0. Particles in the fluid flow cluster if the damping timescale 7 _1 is comparable 
to a timescale characterising the velocity field. In the simulations used in this paper, we 
used a random vector field with a very small correlation time, using the same definitions 
as in [13], where the importance of inertial effects is characterised by a dimensionless 
parameter, which was referred to as e in that work, but which is denoted by rj in 
this paper. It is defined in terms of the correlation function of the velocity gradient 
experienced by a particle with trajectory r(t): 

(ftr) ««>,« (£)<m>). w 

In figure [TJ we showed a realisation of the long-time dynamics. The velocity is 
periodic on the unit square, and was generated from a random stream function 
if)(r,t). The statistics of the stream function are (i[)(r,t)) = 0, (ip{r, t)ip(0, 0)) = 
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Figure 6. The mean number of points in an optimal cover, (Af), as a function of e, 
for various values of a. The slope increases monotonically with a, which ranges from 
a = 0.2 to a = 1 in increments of 0.2. These data, for the random flow model with 
rj = 0.1, show excellent fits to a power law over a wide range of e. 



A 2 exp(— |r| 2 /2£ 2 ) exp(— \t\/r), with £ = 0.1, r small and A chosen such that rj = 0.1. 
The correlation dimension for this value of rj is D 2 ~ 1.76 [H]. 

We examined whether the mean value of the optimal covering number, (A/"), shows 
a power-law dependence upon e. The data for rj = 0.1, shown in figure El are a good fit 
to a power- law for values of a as low as 0.2. At very small values of a, the area of the 
ellipses decreases very rapidly as e — > 0, so that the values of (Af) become too small to 
give reliable results. 

For the inertial particles model we found a number of cases where the covering data 
was not well fitted by a power-law in e. This occurs for small values of a and for values 
of rj where the value of D 2 is close to its minimum, which is D 2 ~ 1.35 at T] ~ 0.35. 
Figure [7] illustrates the case where the fit to a power-law was the least good. 

In figure |8] we exhibit the /3(a) curves for this systems at several different values 
of the inertial parameter t]. Fractal attractors of dynamical systems typically have a 
local structure which is a Cartesian product of a Cantor set and a line. Following the 
discussion in section |2] we therefore expect that the exponent /3(a) should be given by 
equation ([6]), that is /3(a) ~ /3_(a). We find, however, that /3_(a) is not a very good 
approximation, and figure |8] shows that different (3(a) curves may be obtained for cases 
with the same fractal dimension (these arise because D 2 has a minimum with respect to 
varying rf) . In section [2] we also suggested that the generalised dimension /3(a) might 
be higher than /3_(a) when a < | because the unstable manifold is curved. Figure El 
does not show evidence of any discontinuous change at a — \. 

We also investigated /3(a) for fractals generated by two other dynamical processes. 
We simulated diffusion limited aggregation (DLA, one realisation of which was 
illustrated in figure [2]), in the usual manner [?]: points make a random walk on a 
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Figure 7. Same as figure[51 except rj = 0.4, where Z?2 ~ 1-36. In this case the covering 
data are not a good fit to a power-law in e for a < 0.75. 




Figure 8. (3(a) for inertial particles in a random velocity field, using a model defined 
in |14) . The curves are labelled according to the parameter r\ which quantifies the 
importance of particle inertia. 



lattice until they contact the cluster, after which they are frozen and become part of the 
growing set. In the case of isotropic diffusion, the correlation dimension of the resulting 
cluster is D 2 ~ 1.71. The values of the slopes (3 extracted from least-squares fits similar 
to those in figure [6] are plotted in figure M as a function of a. 
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Figure 9. /3(a) for a diffusion limited aggregation (DLA) cluster with 2.7 x 10 5 
particles, and for the Sinai map with a = 0.35. 

We also investigated the fractal generated by the Sinai map, defined by 

x n+ i = x n + y n + a cos(27n/ n )mod 1 

Vn+i = %n + 2y„mod 1 

with a = 0.35, which has an attractor with correlation dimension D 2 ~ 1.62 
curve for this map is also shown in figure [9j 

4. Sierpinski substitution fractals 

The examples of dynamical fractals which we considered in section [3] are all multifractal 
sets, and it is desirable to investigate (3(a) for a model which is a simple fractal, avoiding 
the complications that arise when dealing with multifractal sets [15]. Here we construct 
a class of generalisations of the Sierpinski carpet, which are simple fractals rather than 
multifractal measures. The construction that we use is closely related to one proposed 
independently by Bedford [16] and McMullen [17]. We show that different elements 
from this class of sets can have different (3(a) functions despite having precisely the 
same value of D 2 = /3(1). 

We generate an approximation to a fractal set by a hierarchical process consisting 
of n generations. We generate a set of M n points, where M is an integer, as follows. 
The points Xk lie in the unit square [0,1] ® [0,1], and have coordinates of the form 
(xi,yi) = (i/N^j/Ng), where Ni, N 2 are positive integers satisfying NxN 2 > M. 

We define a 'masking matrix' F with elements as an Ni x N 2 matrix which has 
elements which are either 1 or 0, and we let M be the number of non-zero elements of 




(9) 

. The /3(a) 
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Figure 10. Examples of fractal sets generated by the construction defined in section|4] 
In all of these examples, N i = N 2 = 3. The Cantor set is then defined by listing the zero 
elements of the masking matrix (the labels correspond to the cases considered in figure 
HID: (A) F 22 = 0,D = In8/ln3, (B) F n = 0, D = In8/ln3, (D) F 21 = F 22 = F 23 = 0, 
D = In6/ln3, (E) F 13 = F 31 = F 33 = 0, D = In6/ln3 . 



F. We construct the model set by the following recursive construction. Consider the 
'first generation' set of M points x^, labelled by an index = 1 . . . M, where a point is 
placed at ((i — 1)/Ni, (j — 1)/N 2 ) if i^j = 1. At the next generation each of these points 
is replaced by a set of M points, based on a lattice with spacings N{ 2 and A^ 2 in the 
x and y directions respectively, which are selected by the criterion that Fij = 1. In 
general, after n generations each point Xk is replaced by M points with positions Xk>, 
where k' is an index of the (NiN 2 ) n+1 points, with positions 

, (i - 1 j - 1 \ 

= ** + ^^p' J^nTT ) ■ (10) 

A point labelled by added to the set if and only if = 1. 

As an example, consider the case where N\ = N 2 = 3 and where F 22 = is the 
only element of F which is equal to zero, so that M — 8. Iterating this construction 
gives a version of the Sirepinski carpet set, illustrated in figure [TOlA. with dimension 
D = In8/ln3. By varying the zero elements of the masking matrix, we can generate 
many other Cantor sets, some examples of which are illustrated in the other panels of 
figure [TUl The resulting sets are clearly simple fractals (as opposed to multifractals). 
By making other choices of the masking matrix we can construct other Cantor sets with 
dimension 

D = 21 " M . (11, 

This construction allows us to create distinct fractal sets with exactly the same 
dimension, such as in panels A and B or panels D and E of figure [TUl Moreover, by a 
suitable choice of the masking matrix, we can generate fractal sets which are Cartesian 



Anisotropic covering of fractal sets 



12 



products (such as figure [TUP), as well as those which are not (such as figures ITUR. B, 
E). 

We investigated the function (3(a) for sets which are produced by this generalised 
Sierpinski construction. The results are illustrated in figure [TTJ for seventeen sets 
produced using a 3 x 3 masking matrix. The key at the right hand side of the figure 
indicates the pattern of deletions in the masking matrix, ordered by the number of 
deleted points. 

First we discuss those sets which are a Cartesian product. These include examples 
D, I, L, P and Q in figure [TTJ The simplest example is set D in figures ITUl and ITT1 For this 
set, equation fl6]) predicts that (3(a) = D — 1 +a, which shows excellent agreement with 
figure [TTJ By setting Ni — N 2 — 3 and F 12 = F 2i = F 23 = F 32 = F 22 = we produce 
a set with dimension D = In4/ln3 = 21n2/hi3, which is a Cartesian product of two 
Cantor sets of dimension D x = D y = D/2 = In2/ln3. This is example I in figure [TTJ 
Example L is closely related: this set is similar to example I, rotated by 7r/4. These data 
show quite poor agreement with the prediction from equation (J5l), from which we expect 
(3(a) = D(l + a)/2 (but good agreement with each other). The other two examples 
in figure [TT] which are Cartesian products are very simple: P is a Cartesian product of 
a line and a point, and Q is the product of a Cantor set of dimension In2/ln3 and a 
point. For these examples there is excellent agreement with the predictions of equations 
(EJ) and ([5]), which indicate straight lines of slope unity and In2/ln3 respectively. 

Figure [TT] also shows (3(a) in cases where the set is not a Cartesian product. These 
differ from the data for the Cartesian product sets. They can be organised into sets 
which have apparently identical (3(a) curves. In most of the cases examined in figure 
[TT1 sets which have the same value of M (and hence of D 2 ) have (3(a) curves which are 
identical, to within numerical fluctuations. Examples of such groups are (M = 8 : A, B), 
(M = 5 : G, F), (M = 4 : J, K), (M = 3 : N, O). Note, however, that for M = 4, set M 
is not a Cartesian product and yet has a (3(a) curve which is clearly different from sets 
J and K. 

We noted that examples I and L in figure [HI which are Cartesian products of 
two one-dimensional Cantor sets, had (3(a) functions which show quite poor agreement 
with the expected result, equation ([5]). It appears possible that this anomaly might 
arise because these sets are a degenerate case, where D x = D y . We therefore also 
considered two examples which are a Cartesian product of two Cantor sets with different 
dimensions, namely N\ = 4, N 2 = 3, with non-zero elements Fn = F13 = F41 = F43 = 1 
(which is the set illustrated in figure HJ), and Ni = 5, N 2 = 3, with non-zero elements 
Fn = F13 = F51 = F53. These are Cartesian products of two one-dimensional 
Cantor sets with dimensions D x = | (4 x 3 case) or D x = In2/hi5 (5x3 case), 
and D y = In2/hi3. In figure [121 we show the dependence of (Af) upon e on a double- 
logarithmic scale, for the set illustrated in figure HJ For each value of a, including a = 1, 
the plots appear to have an oscillation about a straight line, which makes it difficult to 
determine accurate values for (3(a). In figure (TSlwe show our best estimates for (3(a) for 
these three fractal sets, compared with the theoretical prediction in the 4x3 case, given 



Anisotropic covering of fractal sets 



13 




Figure 11. f3(a) for the generalised Sierpinski model. The curves are labelled by a 
key giving the zero elements of the masking matrix in white. 



by equation (J5J). There are substantial deviations from the theoretical curve, but note 
that these are no larger than the errors in determining the fractal dimension D 2 = (3(1) 
for these sets. We conclude that although the agreement with (jSJ) is poor, this is due to 
difficulties in fitting the data, and there is no persuasive evidence that ([5} is incorrect. 

Section [2] concluded with an argument which suggests that when the set is locally 
a Cartesian product of a line and a Cantor set, the function /3(a) may have a different 
behaviour when a < \ if the lines are curved. The data on the clustering of inertial 
particles (reported in section [3]) were not a good fit to equation (jBJ), and they did not 
show a clear signature that a = | is a critical point. A more controlled test was 
made with the Sierpinski fractals considered in this section. The sets were modified by 
applying a shift to all of the y-coordinates: 

x = (x, y) — > x' = (a/, y') = (x,y + C sin(27rx)) . (12) 

The sets which are reported upon in figure [TT1 were deformed according to this rule, with 
C — |. The (3(a) exponents remained unchanged for all of these sets (within numerical 
fluctuations comparable to those in figure [TT]) . 
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Figure 12. For the set illustrated in figure 0] it is difficult to fit the exponent /3(a) 
because of an oscillation in the dependence of In (TV) upon lne. 
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Figure 13. /3(a) for three cases of the generalised Sierpinski model which are Cartesian 
products of Cantor sets, based upon 3 x 3, 4 x 3 and 5x3 deletion matrices. 



5. Concluding remarks 

This paper has reported the first systematic study of fractals using a set of covers 
which become more anisotropic as they are made smaller. The growth of the anisotropy 
is described by a parameter a, and we characterised the efficiency of covering by a 
generalised dimension f3(a). We found that different sets with the same correlation 
dimension D 2 can have different (3(a), as was demonstrated very clearly for the simple 
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Sierpinski substitution fractals which we considered in section HI 

In the case where the fractal is locally a Cartesian product of a line and a Cantor 
set, a heuristic argument (similar to that given in [8] and subsequent papers) suggest 
that (3(a) should be given by We find that this expression works well for simple 
model sets of the type considered in section 01 and for the Sinai map, where the attractor 
is locally a Cartesian product of a line and a Cantor set. In the case of clustering of 
inertial particles, however, we find that equation fl6]) does not give a good approximation 
to (3(a). In cases where the Cantor set does not have a Cartesian product structure, 
we were not able to derive an expression for (3(a), and we find persuasive evidence that 
this function is non-universal. 

In three or more dimensions there may also be a tendency for particles to accumulate 
on filamentary structures, as well as on planes. This could be characterised by defining 
the exponent (3 as a function of two parameters, ati, «2 defining a covering ellipsoid with 
principal axes e, e ai and e" 2 . 

This study was partly motivated by the desire to understand the inhomogeneity 
of light scattering from fractal distributions of particles, but our investigations indicate 
that this is a very complex problem, which will be considered in a separate publication. 
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